Interacting bosons in a nearly-resonant potential well 
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We establish that the ability of a localized trapping potential to bind weakly-interacting bosons 
is dramatically enhanced in the vicinity of the threshold of formation of the single-particle bound- 
state of the trap. Specifically, for repulsive particles and a super-threshold trapping potential the 
equilibrium number of bound bosons and the size of the ground state diverge upon approaching 
the single-particle threshold from above. For attractive interactions and a sub-threshold trap a 
collective bound state always forms for a sufficiently large number of bosons despite the inability of 
interparticle attraction alone to form a two-body bound state. 
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I. RESONANT BINDING 

It is empirically known that a surprisingly large num- 
ber of short-range interactions operating in Nature have 
resonant character: the two-particle attraction is either 
barely sufficient to form a two-body bound state or misses 
the binding threshold by a small amount. In nuclear 
physics the prominent examples include the neutron- 
nucleon and neutron-alpha particle attractions which are 
nearly binding [1]. In atomic physics the same is true for 
interactions between the atoms of ■^TJe J2] and of those 
of the spin-polarized Hydrogen family [3]. At the same 
time two '^He atoms form a weakly-bound dimer whose 
binding energy is four orders of magnitude smaller than 
the depth of the '^He-'^He potential well 

The proximity to the two-body binding threshold can 
be quantified by the s-wave scattering length whose mag- 
nitude significantly exceeds the range of interparticle 
forces, thus implying an effectively long-range interaction 
[3| . This regime is fundamentaily interesting because sys- 
tems consisting of resonantly-interacting particles may 
exhibit a series of effects which are independent of the 
microscopic details of the interaction between the con- 
stituents. One such phenomenon is the well-understood 
Efimov effect [H , which occurs in a system of three bosons 
where resonant two-body forces trigger a formation of 
an arbitrarily large number of loosely bound levels in a 
three-particle system. Conclusive experimental observa- 
tion of the Efimov physics became possible only recently 
in an ultra-cold gas of Cs atoms ,6] because, in contrast 
to other nearly-resonant systems found in Nature, the 
two-body scattering length and thus interparticle inter- 
actions in a cold gas can be precisely controlled by using 
Feshbach resonances Q- This allowed the measurement 
of the dependence of the three-body recombination rate 
on the two-body scattering length which contains signa- 
tures specific to the Efimov effect [8] . 

The goal of this paper is to establish the existence of 
a novel resonant effect which may also be observed in 
a cold bosonic gas: we will demonstrate that the num- 
ber of particles bound by a nearly-resonant potential well 



increases as the well is made more shallow, diverging ex- 
actly at the single-particle threshold. The fact that such 
an effect might exists is implied by variational analysis 
given in our previous work [9|. Similar to the Efimov ef- 
fect the ultimate origin of this counterintuitive behavior 
is an effectively large attraction range of the potential 
well. Whereas the Efimov effect requires an attractive 
part of the physical potential to bind three particles, the 
phenomenon discussed below takes place in the presence 
of interparticle repulsion. 

In the laboratory attractive wells have been realized 
using the optical dipole force [l^, in which atoms are 
attracted to the intense region of a focused laser beam. 
Strong confinement in three dimensions can be obtained 
at the intersection of two beams, each typically having 
diameter of about 15 /Ltm [ll| . The binding properties 
of such traps can be tuned by adjusting the laser power. 
Even tighter traps can be achieved by using holographic 
techniques [l^l, optical superlattices [ij], near- field M, 
and white- light techniques (l6j . 

Additionally, tuning the parameters of such localized 
potential placed next to the classical edge of a Bose- 
Einstein condensate, might allow manipulation of a well- 
defined number of particles [oj; it also provides a new 
context for studying the quantum many-body problem, 
as this involves the states of an interacting boson system. 
The experimental feasibility of such a setup was discussed 
in Ref. [§] as well as reasons why tighly focused potential 
is needed. 



II. THE MODEL 

We consider n interacting bosons of mass m in the pres- 
ence of a well- localized attractive potential U{t). We as- 
sume a system of sufficiently low density so that the range 
of interactions between particles need not be included as 
a parameter. Therefore we study the Gross-Pitaevskii 
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(GP) energy functional [T3| 



' ' 2m 



(1) 



subject to the normalization condition 

ip^d^x = 1 (2) 



Here we have assumed that the many-body wavefunction 
can be written as a product of the single-particle func- 
tions '!/'('')• The parameter g represents the short-ranged 
interaction between the particles. The external potential 
U{r) will be chosen in the form of an attractive "shell" 
of radius a: 



U{r) = -aV6{r - a) 



(3) 



because this is a simple form that will bind particles with- 
out having singular behavior in the wavefunction at the 
origin. The case of more realistic potentials will be dis- 
cussed in Section VI. For the case of one particle, there 
is a single bound state for V > Vc = fi.^/2ma^, and this 
condition will continue to be relevant for the many-body 
problem. 



III. RESULTS FROM AN APPROXIMATE 
WAVEFUNCTION 

The ground-state wavefunction minimizes the energy 
functional We can deduce its general form from in- 
spection. If the ground state is to be bound, it will have 
to be relatively large in the well-localized region where 
U{r) is nonzero. At sufficiently large distances it is de- 
creasing more rapidly than 1/r, to ensure normalizability. 
It will pass from large to small in a way that keeps the 
gradient term small; this suggests "0 ~ 1/r at interme- 
diate distances. The effect of the interparticle repulsion 
(the last term in Eq. ([1])) will be to suppress large values 
of tp. 

In a previous paper [9| these considerations led us to 
study the two-parameter trial wavefunction such that 
■0 = C for r < a, and V'(^) = Caexp(— ar -I- aa)/r for 
larger r. We will briefly recapitulate the results. The 
form of the wavefunction for larger r is the exact wave- 
function for a bound state of noninteracting particles, 
and will remain a good approximation to the interact- 
ing wavefunction at large r, where V'^ will be negligibly 
small. The wavefunction extends to a distance of order 
1/a, which will be large for large n and also for V close 
to Vc- Thus we are particularly interested in the case of 
small a. Integration of Eq. Q gives 



(4) 



We note that for small a this reduces to « a/a^, 
rather than (i.e. the reciprocal of the "volume" of the 



wavefunction). This implies that the two length scales 
a and 1/a both play roles; the limit a cannot be 
taken. Eq. ^ leads to 

-2g{n - l)aae'^°"'Ei (4aa)C"* (5) 

where Ei (x) is the exponential integral ^1§\ . 

Using Eq. (|4]) to eliminate mention of a, the leading 
terms of Eq. ^ are 

^ = iV.-V)C^+.ia^V. + lgin-l))C^ (6) 

Minimizing this expression with respect to C gives an es- 
timate for the energy. Assuming g > (i.e. the particles 
repel each other), gives for V >Vc 



E : 



n{V - Vcf 



K + 2.g(n- l)/37ra3 



(7) 



For large n this becomes independent of n, and the cor- 
responding value of a is small. The energy at n ^ oo 
is greater or less than the energy at n = 1, depending 
on whether g is greater or less than g^ = 3Tra^Vc/2. For 
g less than gc, this would seem to indicate that the po- 
tential can bind an infinite number of particles (with a 
very swollen wavefunction); however, including the last 
term of (O gives a shallow minimum for g < gc, and then 
for large n the state is unstable against the unbinding of 
particles (i.e. it is within the continuum of a state with 
smaller n). Our main goal is to establish whether this 
minimum is real, and if this is the case, what is its be- 
havior as the single-particle threshold is approached from 
above, V ^Vc + 0. 

The two regimes for the interparticle strength (the 
cases g > gc and g < gc) are in essence a small-n ef- 
fect. The ground-state energy for two particles is greater 
or lesser than the energy for one, depending on whether 
the interaction parameter is large or small. 

This shows that for V > Vc there is a rt-body bound 
state; however, the total energy for large n becomes inde- 
pendent of n, so that the chemical potential decreases to 
zero. This occurs because for small a the single-particle 
wavefunctions extend to large distances from the attrac- 
tive center, making it correspondingly unlikely that two 
particles will be at the same place, so that the interpar- 
ticle interactions can be effectively made as small as one 
likes (by tuning a). The effect depends on the way the 
size of the ground state scales with the number of bound 
bosons, and does not occur in lower dimensions ^9|]. 



IV. GROSS-PITAEVSKII EQUATION 

The approach just outlined has the advantage that the 
mathematics is transparent, and yields results for the 
energy and wavefunction in which the role of the vari- 
ous parameters in the problem is explicit. It is a robust 
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Figure 1: Dependence of the ground-state energy on the num- 
ber of particles n, for super-threshold well, V — 1.5, and 
various values of the interaction parameter g = 2,4, 6, ...20 
(bottom to top). The solid line is the numerical solution to 
Eq. ([SI substituted into Eq. Q; the dotted line is the mini- 
mum value of Eq. ([5]). Dimensionless variables defined in the 
main text are used. 

method that gives an upper bound on the exact ground- 
state energy. However, the results are only as good as the 
trial wavefunction. Therefore we will find the true mini- 
mum of Eq. ^ by numerical means. This gives the best 
possible wavefunction that can be written as a simple 
product of single particle wavefunctions. In the process 
it verifies that the simple variational function is qualita- 
tively correct, and puts the theory of resonant binding on 
solid footing. Additionally we will show that non-trivial 
cooperative bound states also form for a sub-critical well 
in the presence of a weak interparticle attraction. 

In what follows we will choose a = 1 and Vc — 
h^/2ma^ = 1, which sets the length and energy scales 
for the problem. 

Taking a variational derivative of Eq.([T]) gives the GP 
equation [l7j 

1 

TT (W^) - - 1) + 7^^ = eV- (8) 

where 7 = g{n — 1). Unlike the case of the linear 



Schrodinger equation, the eigenvalue e for which the GP 
equation ([5]) has a normalizable solution is not directly 
proportional to the minimum value of Eq. ([1]); however, 
they are related by 

E = ne-^gn{n^l)B, (9) 

cso 

B = J il^'^inr'^dr (10) 


We solved Eq. ([8]) by integrating from r = with ini- 
tial conditions ip = ^, dip/dr — and chosen values for V, 
7, and e. The integration was in the form of a power-law 
representation of ip{r) for r < 1, which avoided the prob- 
lems coming from the singular form of the GP equation at 
small r and gave accurate values for ip and its derivative 
at r = 1; thereafter, we used second-order Runge-Kutta 
integration with equal-sized steps in the variable y — \[r. 
In general this gave a function that diverged to plus or 
minus infinity at large r. By tuning e we could find a 
case for which ip remained small to chosen large values 
of r. This condition defines the eigenvalue e to high ac- 
curacy. Once a normalizable solution had been found, it 
was rescaled to satisfy the normalization condition 
This also changes the value of 7, but not V or e; only at 
the end of the calculation do we find what value of the 
particle number n the solution corresponds to. 

Figure 1 shows how the ground-state energy varies with 
the number of bosons n, and compares the numerical so- 
lution of the GP equation, Eqs. ©-([H]), to the variational 
result (the minimum value of Eq. (O). The parameter 
V = 1.5, and g = 2, 4, 6, 20. The GP solution is al- 
ways below the variational prediction, because it uses a 
better wavefunction. 

Although the variational results differ quantitatively 
from the solution to the GP theory, qualitatively they 
are similar. We offer the following observations: 

*The energy is negative for all n, so that the attractive 
potential has a bound state with an infinite number of 
particles. 

*For g < gc, the energy has a very shallow minimum as 
a function of the number of particles. This implies that 
it will be difficult to attach a definite number of particles 
to a well by placing it in tunnelling contact with a source 
(such as a condensate droplet): there is no equilibrium 
at all, or the number of particles bound will be large and 
very susceptible to small perturbations. The minimum in 
E{n) is even less pronounced in the GP solution than it is 
in the variational result. We have indicated the apparent 
position of the minima, using the graph marker. 

*For g > (^cj the energy increases with the number 
of particles. Then it is possible to bind a well-defined 
number of particles by letting the system equilibrate with 
a source. However, this scheme will only work well for 
a small number of particles: the plateau at large n will 
make the system insensitive to the number of particles. 
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Figure 2: How rix, the value of particle number n that min- 
imizes the energy, varies with the strength V of the binding 
potential. The interaction parameter was taken to be <; = 2. 
The dotted line corresponds to the variational result. 

In the variational study, we found a counterintuitive 
behavior for the value at which the minimum of the 
total energy occurs: when g is below the critical value gc, 
Ux becomes large, diverging as V approaches the binding 
threshold. This means that an arbitrarily large number 
of particles can be stably bound by making the potential 
weak! We have verified that this behavior is also exhib- 
ited by the solution to the GP theory. Here it is useful 
to observe that we are trying to find the minimum of the 
total energy both with respect to the parameter n and 
the normalized wavcfunction 7p. Imposing the condition 
that the derivative of Eq.(l) with respect to n (with ij} 
fixed) is zero, and combining the outcome with Eq. ([5]), 
we obtain a relationship 

g^2\e\/B. (11) 

Imposing this condition proved to be easier than trying to 
minimize E over n. We constructed Figure 2 by finding 
normalized solutions to the GP equation ([5]) for various 
values of 7, and iterating until Eq. pT|) holds. 

Although our work demonstrates the existence of the 
effect of resonant binding only within the GP approx- 
imation, we expect that it is also present in the exact 
theory. As a nearly-resonant well is made more shal- 
low approaching the single-particle binding threshold, the 
effective range of the central attraction increases, and 
the well binds an increasingly larger number of bosons. 
Therefore the microscopic details of the well are not ex- 
pected to be of qualitative significance. The large range 
is responsible for the survival of the effect in the pres- 
ence of weak short-range repulsion between the bosons 



because it allows the particles to take advantage of the 
central attraction while minimizing mutual contact in- 
teractions [Toj . The interactions become more important 
as the well is made deeper, thus explaining the existence 
of the minimum in Fig. 2. For a sufficiently deep well 
the number of bound particles Ux increases as the well 
deepens which is the normally expected behavior. 

Additionally, we argue that our results are immune to 
the effects of correlations neglected in the GP approxi- 
mation. Indeed, the GP theory can be viewed as a varia- 
tional method. Therefore the exact energy will be lower 
than our calculations indicate. What they do show is 
that near resonance the wavcfunction can accommodate 
an unexpectedly large number of particles, because each 
particle manages to sample the attractive well without 
interacting with the other particles too strongly. The 
effect of correlations would be that the particles evade 
each other more effectively, and thus enhance the effect 
we predict. 

The issue of validity of the GP theory was also ad- 
dressed in Ref.0 where the same problem was studied 
in general dimension. Specifically, in one dimension a 
version of the problem in question was solved exactly 
by using the Bethe ansatz methods and the results were 
compared with the GP solution. Although the exact 
ground-state wavefunction is indeed not of the product 
(Hartree) type, in the limit of large particle number the 
exact ground-state energy turns out to be nearly identi- 
cal to its GP counterpart. This example represents "the 
worst case scenario" as the effect of correlations is ex- 
pected to be strongest in one dimension. Therefore we 
conclude that it is very likely that the GP theory is fully 
adequate in three spatial dimensions. 



V. WEAKLY ATTRACTING PARTICLES 

The variational analysis also predicts an interesting 
phenomenon for the case of weakly-attracting particles 
in the presence of a potential that is too weak to bind a 
single particle (g < and ^ < 1). The second term of 
Eq. ^ is negative for sufficiently large rt, pointing to the 
existence of a many-body bound state. For the particular 
case V = 1, the critical value of the number of particles 
needed is given by ric = 3Tr/2\g\. When V < 1, the min- 
imum energy can be positive (indicating a metastable 
state), but becomes negative for sufficiently large n. The 
possibility of a cooperatively bound state of this sort was 
first pointed out by Migdal [2^. 

Figure 3 shows the results of numerical solution of the 
GP theory for the case of weakly attracting particles in 
the presence of a y < 1 potential. In agreement with the 
variational argument, collective bound states are seen to 
form for a sufficiently large number of bosons. Specifi- 
cally, the farther away is the binding threshold, the larger 
number of bosons is required to form a bound state. 
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Figure 3: The ground-state energy of weakly-attractive 
bosons as a function of their number n, for various values 
of the binding potential amplitude V. The curves are all for 
the case g — —0.1, and the successive curves from left to right, 
are for V = 1.00, 0.98, 0.96, 0.94, 0.92, 0.90. 

VI. BEYOND THE MODEL 

The delta-shell potential cannot be achieved experi- 
mentally. Then one may ask how our results are relevant 
to the real world. 

The variational results are quite general for a local- 
ized external potential. Since the trial wavefunction is 
constant for r < a, the external potential only enters 
through its average value 
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The solution to the GP equation will be affected slightly. 
We can anticipate the effects by comparing the delta shell 
problem to that of the attractive square well U{r) = —W 
for r < a. Outside the well the GP equation is the same 



for the square well and the delta shell, and so the solu- 
tions (for given 7 and e) are the same; the only question 
is how the value for W corresponds to that for V. The 
condition that relates them is that the values of the log- 
arithmic derivative of the wavefunctions should be the 
same at r = a. For the delta shell, the solution to Eq.® 
(for 7 = 0) is V' = sinh(-\/|ejr)/r, and is very nearly 
constant for the cases of interest; however for the square 
well the solution inside (again for noninteracting parti- 
cles) has the form ip — C sm{^yW — \e\r)/r. The loga- 
rithmic derivative of this function takes on all values and 
there is an infinite set of values of W for which the match 
can be achieved. The effect of the particle interactions 
will be to somewhat suppress the accumulation of parti- 
cles (places where ip is large) and to push the consecutive 
bound states to larger values of W. 

However, the feature of interest in the resonantly 
bound system is the number of particles that can be 
bound, and this is a property of the wavefunction at large 
distances. Thus we can expect that the the differences be- 
tween different external potentials can be absorbed into 
the definitions of Vc and g. 



VII. CONCLUSIONS 

To summarize, our work demonstrates that a super- 
threshold localized trapping potential can bind an unex- 
pectedly large number of weakly-repulsive bosons while 
a sub-threshold trap can trigger a formation of a cooper- 
ative bound state of weakly-attractive bosons even when 
two-body attraction is insufficient to form a two-body 
state. These effects are insensitive to the microscopic de- 
tails of the trapping potential and have their origin in the 
large effective range of nearly-resonant trapping poten- 
tial. We hope that these results will be experimentally 
tested in the near future. 
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